Generation of Werner-like stationary states of two qubits in a thermal reservoir 
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The dynamics of entanglement between two-level atoms immersed in the common photon reservoir at finite 
temperature is investigated. It is shown that in the regime of strong correlations there are nontrivial asymptotic 
states which can be interpreted in terms of thermal generalization of Werner states. 
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I. INTRODUCTION 

Entanglement generation by the indirect interaction be- 
tween otherwise decoupled systems has been discussed in the 
literature mainly in the case of two-level atoms interacting 
with the common vacuum. The idea that dissipation can cre- 
ate rather then destroy entanglement in some systems, was put 
forward in several publications 01 HI HI 0] • In particular, the 
effect of spontaneous emission on destruction and production 
of entanglement was discussed JH 0, 0- When the two 
atoms are separated by a distance small compared to the ra- 
diation wave length, the collective properties of two-atomic 
systems can alter the decay process compared with the sin- 
gle atom. It was already shown by Dicke H that there are 
states with enhanced emission rates (superradiant states) and 
such that the emission rate is reduced (subradiant states). (In 
the case of multi-atomic systems, similar collective effects de- 
scribed by Tavis-Cummings interaction [10] can generate col- 
lective multiqubit entanglement |1 1]). When the emission rate 
is reduced, two-atom system can decohere slower compared 
with individual atoms and some amount of initial entangle- 
ment can be preserved or even created by the indirect inter- 
action between atoms. The analogous effect of production of 
entanglement was studied in a system of two-level atoms in- 
teracting with a squeezed vacuum. In that case, the squeezed 
vacuum is a source of nonclassical correlations which are es- 
sential for creation of entanglement [ 12]. The general condi- 
tion under which entanglement can be induced by interaction 
with environment was also discovered 

The case of two atoms immersed in a common thermal 
reservoir was also investigated ffl fl5i [l6ll . In particular, Be- 
natti and Floreanini [ 14] have discussed the dynamics of two 
independent atoms interacting with the reservoir of scalar par- 
ticles at finite temperature and found the interesting behavior 
of the system. If the atoms are at finite separation, there is a 
temperature of the reservoir below which entanglement gen- 
eration occurs. Moreover, for the vanishing separation, the 
entanglement thus generated persists in the asymptotic state. 

In the present paper, we study the similar model but from 
a different perspective. We consider the system of two-level 
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atoms interacting with the photon reservoir at fixed tempera- 
ture. As in the vacuum case, the collective properties of the 
atomic system can alter the decay process compared with the 
single atom. There are states with enhanced emission rates 
and such that the emission rate is reduced. The important 
example of the latter is the singlet state \a) i.e. antisymmet- 
ric superposition constructed from energy levels of considered 
atoms. As we show, in the regime of strong correlations the 
singlet state is decoupled from the environment and therefore 
is stable. In our considerations we focus on this case and in- 
vestigate the stationary asymptotic states of the system. In 
turns out that the asymptotic states are parametrized by the fi- 
delity F of the initial state p with respect to the state \a) (or 
the overlap of p with singlet |a)) and the temperature T of 
the photon reservoir. We identify the asymptotic states as the 
thermal generalization of Werner states i.e. mixture of singlet 
state and Gibbs equilibrium state at the temperature T with 
the appropriate probability. (The standard Werner states lfl7ll 
are recovered in the limit of infinite temperature). Depending 
on the initial fidelity some of the asymptotic states are entan- 
gled. We calculate the amount of the asymptotic entanglement 
using a concurrence as its measure. We also show that if the 
initial fidelity is greater or equal to 1/2, then for every finite 
temperature of the reservoir the asymptotic entanglement is 
non-zero. On the other hand, for fidelity less then 1 /2 there 
is the critical temperature below which the asymptotic states 
are entangled whereas for higher temperatures all stationary 
states are separable. 

Since the initial states with the same fidelity can be separa- 
ble or entangled, the dynamics given by the interaction with 
the reservoir manifests differently with respect to the entan- 
glement properties of the system. It can create entanglement 
when the initial states are separable, disentangle initially en- 
tangled states, preserve some part of initial entanglement or 
even increase the initial entanglement of some states. We 
show on specific examples that this behavior of the dynamics 
really occurs. In particular, pure product states with orthog- 
onal factor vectors become entangled in every temperature, 
whereas for other product states there is the critical temper- 
ature for creation of entanglement. Similar phenomenon can 
be observed if the initial state is mixed. For example in the 
case of Gibbs equilibrium state at the temperature 7b which 
differs from the temperature T of the reservoir, the creation 
of asymptotic entanglement happens if To is much higher then 
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T. On the other hand, entangled initial states with maximal 
amount of entanglement can preserve some initial entangle- 
ment or can disentangle completely. It is worth to stress that 
this different behavior with respect to the thermal noise can 
happen for locally equivalent initial states, so by performing 
only local operations one can protect as much entanglement as 
possible. When the initial states are not maximally entangled, 
then the thermal noise can in some cases create entanglement 
which adds to the initial one, although during the evolution 
the purity is decreasing. 



II. TWO QUBITS DYNAMICS 

We start with the sketch of the derivation of the dynami- 
cal equation describing the evolution of the system of atoms 
immersed in the thermal reservoir (for details see e.g. B18ID . 

Consider two-level atoms A and B with ground states \0)j 
and excited states \l)j (j — A,B), interacting with photon 
reservoir at temperature T. The dynamics of the combined 
system consisting of the atoms and the quantum electromag- 
netic field is given by the Hamiltonian 

H SF =H +H I 

where Hq is the sum of free Hamiltonians of the atoms and 
the field and Hj describes the interaction between the atoms 
and photons in the electric dipole approximation. Since we 
are interested in the dynamics of the system of atoms, we take 
the partial trace over field variables and find that the reduced 
density matrix of the atomic system satisfies some integro- 
differential equation which can be simplified by employing 
the Born approximation (the interaction between the atoms 
and the field is so weak that there is no back reaction effect 
of the atoms on the field). Under this approximation, the evo- 
lution of the density matrix depends on the first- and second- 
order correlation functions of the field operators in the ther- 
mal equilibrium state at temperature T. Employing further the 
rotating-wave approximation in which we ignore all terms os- 
cillating at higher frequencies and assuming that the correla- 
tion time of the photon reservoir is short (the Markov approx- 
imation), we arrive at the result that the influence of reservoir 
on the system of atoms can be described by dynamical semi- 
group [19] with Lindblad generator L = —i[H, ■]+ Lq, where 



H = J E °3 + I %*<0* 

~ j=A,B j,k=AM 



(II. 1) 



and 



l dP=2 L r l jk (2oLpo k + -o k + oLp-po k + o J 



■ j,k=A,B 



(II.2) 



+ rj, (2a j + pa k _ a V +P - paV; 



Here 



oj_ = o±®l, of = l®a±, of =o 3 (8)1, af = ll<g)a3 



In the Hamiltonian (III, lb . co is the frequency of the transition 
|0)y — > (j = A, B) and Q.ab — &ba — Q describes inter- 
atomic coupling by the dipole-dipole interaction. On the other 
hand, dissipative dynamics is given by the generator ill.21 
with 



where 



Ijk (! + (»», 



-p«> 



-A . 



1 



,-Pco' 



is the mean number of photons, and 



Jaa = Jbb = Yo, Jab = Yba = Y 



(H.3) 



In the above equalities, yo is the single atom spontaneous 
emission rate, and y = Gyo is the collective damping constant. 
In the model considered, G is the function of the interatomic 
distance R, and G is small for large separation of atoms. On 
the other hand, G — > 1 when R is small. 
The master equation 



d P r„ 



(II.4) 



giving the time evolution of a density matrix of the system 
of two-level atoms can be used to obtain the equations for its 
matrix elements with respect to some basis. To simplify the 
calculations one can work in the basis of collective states in 
the Hilbert space C 2 <8 C 2 [18], given by product vectors 



\e) = \l) A ®\l) B , |s) = |0>A<8|0)i 
symmetric superposition 

k) = -^(|0>A®|l>H + |l>A®|0>lO 

and antisymmetric superposition 
1 



V2 



(|1)a®|0) b -|0)a®|1)b) 



(H.5) 



(H.6) 



(H.7) 



In the basis of collective states, two-atom system can be 
treated as a single four-level system with ground state ex- 
cited state \e) and two intermediate states \s) and \a). From 
( 1H.41 i it follows that the matrix elements of the state p with re- 
spect to the basis |e), \s), \a), \g) satisfy the equations which 
can be grouped into decoupled systems of differential equa- 
tions. So for diagonal matrix elements we obtain 



dp e 



dt 

dPss 
dt 

dpaq 

dt 

dpgg 
dt 



(YO - Y) («) Paa - 2y ( 1 + («) ) Pee + (YO + Y) («) P.v.v 

-(Yo+Y)[(l + (n))Pee- (1 +2<n»p„+<B>p„] 
-(Yo- Y)[(l +2{n))p aa - (1 + (n))p ee - (n)p gg } 
(Yo - Y) (1 + (») )Paa + (Yo + i)p ss - 2yo («) Pgg 



(118) 
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On the other hand, the elements p es and p sg are connected by 
the equations 



dPes 

dt 
dps 



{jQ+j)(n)p sg --[(l+2(n))j 

+ (3+4<B))7b-2/((n-n)]p B 
1 



(H.9) 



= (Yo + Y)(l + («»P«- 5KI +2<n))y 
+ (l+4(n))y + 2/(co + n)]p s? 
and similarly, the elements p^ and p ag satisfy 

^p=-( 7b -Y)< B )p <w +I[(l+2< B ))Y 



t/p 



(3 + 4(«))y -2/(co + n)]p 
1 



(11.10) 



dt 



-(Yo-Y)(l + («))pe 



(l+4(ra))Yo 



-(l+2(«))y+2/(co-ii)]p ag 



Finally 



and 



t/p e 



dt 



dpsi 
dt 



= -[(l+2(«))y + 2/co]p 



= -[(l+2(n))7b + 2/a]p« 



(11.11) 



(11.12) 



The equations for the remaining matrix elements can be ob- 
tained by using hermiticity of p. 

From the equations (III. 8b it follows that similarly as in the 
zero temperature case (see e.g. II 1810 . the system of atoms 
prepared in the symmetric state \s) decays with enhanced rate 
Yo +Y> whereas antisymmetric initial state \a) leads to the re- 
duced rate yo — y. In the limiting case of strongly correlated 
atoms we can put y = yo, so the state |a) is completely de- 
coupled from the photon reservoir. One can also check that 
the master equation ( 111.4) describes two types of time evolu- 
tion of the atomic system, depending on the relation between 
y and yo. When y < yo, there is a unique asymptotic state of 
the system, which is the Gibbs state 



Pp 



/tr 



2 H ~ 3 

j=A,B 



(11.13) 



The state ( III. 13b is separable and describes thermal equilib- 
rium of atoms interacting with photon reservoir. In the regime 
of strong correlations, y = yo and we show that there are non- 
trivial asymptotic stationary states which can be parametrized 
by matrix elements p aa of the initial state. 



III. 



STRONGLY CORRELATED QUBITS AND 
NONTRIVIAL ASYMPTOTIC STATES 



When y = yo, equations (IH.8I 1 - ( 111.12) simplify and one can 
check that the solutions of (1H.9) - (III. 1 2b asymptotically van- 
ish, so the only contribution to the stationary states p as comes 
from the matrix elements p aa , p ss ,Pee and p gg . Notice that 



dp ai 
dt 



= 0, SO Paa(t) = Paa(0) = F 



where 



F = (a\p\a) 



is the fidelity of the initial state p with respect to the singlet 
state \a). Hence (p as )«a = F and after a long elementary cal- 
culation, we obtain that 

(Pa.)« = -7ov- (!-*") 



Il(P) 



(III.l) 



(Pas)gg 



1 



(1-F) 



gS ' «(P) 
where 

In the canonical basis 

|1}a®|1>* |l>A®|0)n, |0)a®|1)b, |0) A ®|0) fi 
the non-zero matrix elements of the asymptotic state reads 



(Pas)l 



«(P) 



(l-F) 



(Pas)22 ^MP) (1 - F)+ 2 



(Pas)23 



(Pas) 



44 



g-fr" F 
(l-F) 



(in.2) 



«(P) 



and (p as )33 = (Pas)22- 

The asymptotic state p as defined by dIII.2| i exists for any ini- 
tial state and for the fixed temperature of the photon reservoir 
it depends only on the initial fidelity i.e. p as = p as (F). If we 
define the threshold fidelity Fp by 



„-P<B 



F R = 



(l +e -PtO)2 

then one can check that: 

(1) p as (Fp) equals the the Gibbs state pp, 

(2) for F > Fp, p as (F) equals to the thermal generalization of 
the Werner state 



Wp = (l-p)pp + p\a)(a\ 



(1113) 



where the mixing probability p depends on the fidelity of the 
initial state and temperature of the reservoir : 



(1 +e -P»)2 F _ e -Pco 



(III.4) 



(3) for F < Fp, the state p as cannot be expressed as the Werner 

state dnni ) 
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Notice that in the limit P — > (or T — > oo) we obtain the 
standard Werner state 1 17], hence the dynamical generation of 
such states occurs due to the interaction with the environment 
with maximal noise [20]. The entanglement properties of the 
asymptotic state will be discussed in the next section. 



IV. ASYMPTOTIC ENTANGLEMENT 

We start with the characterization of entanglement of the 
thermal Werner state Wp. The simplest way to do this is to use 
Wootters concurrence B21I1 defined for any two-qubit state p 
as 



C(p) = max I 0, v %\ — y^-2 — V^3 — 



(IV. 1) 



where A-i > \q > A-3 > A4 are the eigenvalues of the matrix pp 
with p given by 

p = 02(^02 p 02(^02 

where p denotes complex conjugation of the matrix p. By a 
direct calculation one obtains that for the states 



/pn 

P22 P23 

p 32 P33 

P44/ 



(IV.2) 



the concurrence is given by the simple function 

C(p)=max(0, 2(|p 23 |-VpTTp47)) (IV.3) 
Since the states ( IIII.3I ) are of the form (IIV.21 i, C(Wp) reads 

2e-P m \ 

' p)\ (IV.4) 




+ e -2Pco 



thermal Werner states are entangled, and for p < po, those 
states are separable. Combining this result with the formula 
( IIII.4t . we obtain the concurrence of the asymptotic states (see 
aisoTlIGl) 



C as = max 0, F 



1 +2coshpco 



(1 



(IV.5) 



Observe that for the fixed temperature of the reservoir the 
asymptotic entanglement depends only on the fidelity of the 
initial state. Moreover, this entanglement is non-zero for all 
initial states with fidelity satisfying 



F>F 



4 + 2coshpco 



(IV.6) 



Obviously, all asymptotic states with fidelity less then thresh- 
old value fp are separable. 



We can also consider the interesting problem of tempera- 
ture dependence of the asymptotic entanglement. Notice that 
for every p > and F > 1/2 

(1-F)>0 



1 +2coshpco 



So the interaction of the atomic system with the reservoir at 
any finite temperature brings all initial states with the fidelity 
F greater or equal to 1 /2 into the stationary entangled states. 
On the other hand, when the fidelity is smaller then 1/2, there 
is the critical temperature i.e such temperature T c that if T < 
T c the asymptotic states are entangled whereas for T >T C the 
asymptotic states corresponding to the same initial fidelity are 
separable. Simple calculation shows that 



CO 



= <p(F) 



where for F < 1/2 



q>(F) 



In 



'VV3-8F + 4F 2 + 3-4F N 
IF 



(IV.7) 



(IV.8) 



Since the states with the same fidelity can be separable or 
entangled, we expect that the system will behave differently 
depending on the initial conditions. More precisely, it can 
happen that unentangled atoms become entangled during the 
evolution and initially entangled states disentangle or remain 
entangled. Next we show that the above possibilities actually 
occur. 



A. Separable initial states 



When the initial state is the pure product state 



|V) = |9>®|V), l<P),l¥>eC 2 



then its fidelity is given by 



F = i(i-i<cp| ¥ )i 2 ; 



(IV.9) 



(IV. 10) 



If we denote |(9|v|/)| = a, the concurrence of the asymptotic 
state corresponding to flIV9t can be computed from the for- 
mula 

, N / (l-a 2 )coshpco-l-2a 2 \ 

^(«)= max (o^ i T ^L j (IY11) 

In the case when the vectors cp and \|/ are orthogonal (for ex- 
ample (p is the excited state of the atom A and \|/ is the ground 
state of the atom B), F is maximal and for every finite temper- 
ature the concurrence of the asymptotic state is non-zero. So 
the interaction with the photon reservoir creates the station- 
ary entanglement between initially unentangled atoms. The 
amount of this entanglement is maximal for the zero temper- 
ature and decreases asymptotically to zero when the tempera- 
ture of the reservoir increases to infinity(FIG. 1). Notice also 
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FIG. 1: Ca S (a) as a function of temperature for the state HV.9b with 
a = 0. 



FIG. 3: Critical temperature versus 7j). 




FIG. 2: C as (a) as a function of temperature for the state i ll V. 91 ) with 
a = 1/4. 



that for a = 1 , C as (oc) trivially vanish. In the general case 
when < a < 1, the creation of the stationary entanglement 
is possible only for the temperatures of the reservoir below the 
critical temperature, which in that case is given by the formula 



7c 

CO 



In 



'v / 3V2a 2 + oc 4 + l+2a 2 



-i -l 



1 



(IV. 12) 



Since T c goes to infinity when a goes to zero, the critical tem- 
perature can be high for some initial states, but the created en- 
tanglement is always maximal for zero temperature and van- 
ishes for T = T c (FIG. 2). 

As the example of mixed separable initial state, consider 
the Gibbs state dll. 1 3b at some temperature Tq ^ T. Notice 
that the fidelity F of that state is given by 



F = 



(1 



(IV. 13) 



and the value of F is always smaller then 1/2. So there is 
the critical temperature T c of the reservoir for the creation of 
asymptotic entanglement and T c depends on Tq. Explicit func- 
tion describing T c can be obtained by combining the formulas 



( lIV.7b and dlV. 1 3b . Since the resulting function is complicated, 
we will not reproduce it here and show only its plot (FIG. 3). 
Observe that the critical temperature is always below the tem- 
perature To, thus in order to entangle the atomic system pre- 
pared in the Gibbs state at temperature 7b, the temperature of 
the photon reservoir have to be much smaller then To- On the 
other hand, the amount of the created entanglement can be 
computed from the formula dIV5b . Similarly as in the case 
of pure initial states, asymptotic concurrence is a decreasing 
function of temperature and vanishes at T c (FIG. 4). 




0.3 0.4 T c /oj 

T/CD 

FIG. 4: C as as a function of temperature for the Gibbs state with 
r /eo = 4. 



B. Entangled initial states 

The collective states \s) and \a) are specific instances of 
states with maximal possible entanglement. In the case of 
two qubits the set of maximally entangled states form a three- 
parameter family of states. As was shown in [22], the corre- 
sponding projectors can be parametrized as follows 
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a 



ce 



't>2 



v- 



. a 2 e i{^ 2 ) 

,2 _ 1 „2 



fc 2 g -i(d 1 -i5 2 ) ^2 



-m 



-ce 



2 

2 



a 



J 

(IV. 14) 



where b 2 = 1 — a 2 , c = a\/l — a 2 and 

ae [0,1], r>i,r> 2 S [0,27t] 
All the states (IIV. 14b have concurrence equal to 1, but 

F=i(l-a 2 )(l-cos(t>i-f> 2 )) (IV.15) 
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Notice that the fidelity F can take all values from to 1, 
depending on parameters a and i3- = "&\ — 13-2- In particular 
F > 1/2 inside the set £ on the plane, given by 



1 a 2 

{ < a < — = , arccos ^= 

V2 a - 1 



< f> < 2% — arccos - 



} 



a 2 -l 
(IV. 16) 

On the boundary of £ , F = 1/2 and outside this set, F < 1/2. 
So all initial states (IIV. 14b with (a,t>i - #2) £ £ remain en- 
tangled asymptotically, for any finite temperature of the reser- 
voir. The asymptotic concurrence can be computed from the 
formula 



C as (a,#) 



2(1 -a 2 ) (sin 2 |coshPco- cost}) - (1 +2a 2 ) 

1 +2coshpco 

(IV. 17) 

for (a,t>) G £ . Notice that for all such initial states the value 
of asymptotic concurrence is smaller then 1, except antisym- 
metric collective state \a) which is stable during the evolution. 
As in the previous cases, the asymptotic concurrence is maxi- 
mal for the reservoir at zero temperature, with the value given 
by 



Cmax — (1 



2, ■ 2* 

■a ) sin 2' 



(IV. 18) 



but in the present case, it decreases to the nonzero minimal 
value C m ; n when the temperature of the reservoir increases to 
infinity, and 



C min = - (a 2 + (l-a 2 )cosr>) 



(IV. 19) 



For the parameters (a,r>) lying outside the set £, there is 
the critical temperature T c given by fl!V.7t and ( |IV.8t for fi- 
delity ( IIV. 15b . For T < T c , the interaction with reservoir pre- 
serves some initial entanglement, but for T > T c , all maxi- 
mally entangled states with such parameters a and $ disen- 
tangle asymptotically. Some initial states (1IV.141 > disentangle 
exactly to the equilibrium Gibbs states at the specific temper- 
ature T, depending on a and i3-. This interesting phenomenon 
can happen when 



(1-a 2 ) (1-cosd) = 



,-pco 



(l+e-N) 2 



p = i (IV.20) 

T 



FIG. 5: Asymptotic entanglement as a function of temperature for 
the initial state JIV.23I >. Here C is the initial entanglement. 



This equation can be satisfied for some finite temperature T 
only if the parameters (a, ft) lie outside the curve 



2a 2 -l 



r> = arccos — — — , a G\0,V3/2] 
2 (a L — 1 ) 



(IV.21) 



and then 



ln- 



( 1 - a 2 ) cos f> + y/2a 2 - 1 + 2( 1 - a 2 ) cos r> 



(l-a 2 )(l-cosr>) 

(IV.22) 

On the curve (IIV.21I ). this temperature is infinite. 

When the initial states are non-maximally entangled, the in- 
teraction with the photon reservoir can in some cases increase 
the initial entanglement. To show that this is possible, con- 
sider the class of states 



(IV.23) 
















X 


~2~ 








z 


1 -x 





Vo 








0) 



where (0,1) are related by inequality 

— <x(l -x) 
4 - V ' 



Note that as the explicit examples of states ( IIV.23I I we can 
consider the pure states 

1*1% - cost] |0) A <g>|l)s + sinri 1 1) A ® |0) B (IV.24) 



with T| S (ti/2, 7t). The states (|IV.231 > are entangled with con- 
currence equal to z, and 

2 2 

For every temperature of the reservoir the asymptotic concur- 
rence is non-zero and is given by 



Cas(z) + 



(l-z)(coshpco-l) 
1 +2coshpco 



(IV.25) 
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So it is greater then the initial concurrence and the maxi- 
mal production of the additional entanglement happens at zero 
temperature (FIG. 5). In that case 

Observe also that when T — > oo, C as (z) — > z, so in the reser- 
voir at infinite temperature the initial entanglement is exactly 
preserved 12011 . 

V. CONCLUSIONS 

We have investigated the dynamics of two-level atoms im- 
mersed in the photon reservoir at finite temperature T. In 
the regime of strong correlations between the atoms, there are 
nontrivial stationary asymptotic states which are parametrized 



by the fidelity F i.e. the overlap of the initial state of the 
atoms with the singlet state and the temperature T. For the 
values of F above the threshold fidelity, these states can be 
identified with thermal Werner states, which are natural gen- 
eralizations of standard Werner states. Depending on F and T, 
the asymptotic states can be separable or entangled. Thus the 
dynamics describes the process of creation of entanglement or 
the phenomenon of disentanglement. Concerning generation 
of entanglement it is worth to stress that there exists a criti- 
cal temperature above which the entanglement cannot be cre- 
ated. Besides, even in the case of generation of entanglement, 
the temperature diminishes its production. The maximal value 
of entanglement is obtained for the case of zero temperature. 
On the other hand, in the process of disentanglement some 
part of initial entanglement can be preserved but we can also 
find such entangled initial states that disentangle exactly to the 
Gibbs equilibrium state. 
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